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We have estimated the critical exponent describing the divergence of the localization length at the 
metal-quantum spin Hall insulator transition. The critical exponent for the metal-ordinary insulator 
transition in quantum spin Hall systems is known to be consistent with that of topologically trivial 
symplectic systems. However, the precise estimation of the critical exponent for the metal-quantum 
spin Hall insulator transition proved to be problematic because of the existence, in this case, of edge 
states in the localized phase. We have overcome this difficulty by analyzing the second smallest 
positive Lyapunov exponent instead of the smallest positive Lyapunov exponent. We find a value for 
the critical exponent v = 2.73 ± 0.02 that is consistent with that for topologically trivial symplectic 
systems. 

PACS numbers: 71.30.+h,73.20.Fz 

I. INTRODUCTION 

Onoda et al.— have questioned whether Anderson transitions in topologically non-trivial systems share the same 
critical properties as those in topologically trivial systems. Obuse et al., 2,3 studied the metal-ordinary insulator 
transition in quantum spin Hall (QSH) systems and found that the value of the critical exponent is the same as 
that of topologically trivial symplectic systems. However, Obuse et al. only studied transitions to insulating phases 
without edge states leaving open the possibility that the critical exponent for the metal-QSH insulator transition 
might be different. Indeed, the critical conductance distributions have been found to be sensitive to the number of 
the edge states^ 

Here, we report an estimation of the critical exponent for the divergence of the localization length at the metal-QSH 
insulator transition. We find a value of the critical exponent that is consistent with that in the conventional symplectic 
class. Our result supports the conjecture that the metal-QSH insulator transition belongs to the conventional Wigner- 
Dyson symplectic class. 

A. Second smallest positive Lyapunov exponent 

In order to estimate the critical exponents for the divergence of the localization length, we have performed a finite size 
scaling analysis of numerical data for Lyapunov exponents^— This involves the estimation of the Lyapunov exponents 
for quasi-one-dimensional systems with effectively infinite length and given cross sections. For extrapolation to the 
two dimensional limit that is of interest here, a strip with a cross section L is considered. For ordinary metal-insulator 
transitions, the standard approach is to analyze the scaling of the dimensionless quantity 

r = r 1 = 7l L, (i) 

which is equal to the product of the smallest positive Lyapunov exponent 71 (precisely defined below) and the cross 
section. (This quantity is the inverse of the so-called MacKinnon-Kramer parameter.) In a metallic phase T usually 
increase with L and in a localized phase usually decreases with L. At the transition point between these phases T 
becomes independent of L and a common crossing point of curves with different L is visible on the appropriate graph. 
However, when applying this method to the metal-QSH insulator transition we run into a problem: T increases in 
both the metallic and QSH insulating phases and no common crossing point is seen (Fig. ^ . This problem, which 
makes a very precise scaling analysis difficult, occurs because of the existence of edge states in the QSH insulating 
phase. We have overcome this problem by analyzing the scaling of the quantity 

r 2 - 7 2 i, (2) 

which involves the second smallest positive Lyapunov exponent. This Lyapunov exponent is much less affected by 
these edge states and a common crossing point is recovered (Fig. [3]) . It has been demonstrated^ for the Anderson 
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FIG. 1: Phase diagram of Z2 network model with (a) PBC and (b) RBC. When RBC are imposed, the insulating phase 
without edge states appears for small p, the quantum spin Hall insulating phase appears for large p, and a metallic phase 
appears between the two insulating phases. The arrow indicates the range of the parameter considered in this paper. 



transition in three dimensional systems in the orthogonal symmetry class that the critical exponent obtained from the 
scaling of higher Lyapunov exponents is the same as that obtained from the scaling of the smallest positive exponent. 



II. CALCULATION OF LYAPUNOV EXPONENTS 



A. Model 



To describe the QSH system, we use a Z2 quantum network model.— ~— The Z2 network model has two controlling 
parameters, p, the tunneling probability (related to the chemical potential), and q, the spin- mixing probability (related 
to the spin-orbit interaction strength). For finite q, the Z2 network model exhibits a metallic phase sandwiched between 
two insulating phases. When periodic boundary conditions (PBC) are imposed, both insulating phases correspond to 
the ordinary insulating phase [Fig.QJa)]. On the other hand, when reflecting boundary conditions (RBC) are imposed, 
the insulating phase located for larger p exhibits edge states and becomes the QSH insulating phase [Fig. [2(b)] . Thus, 
the Z2 network model with RBC shows two types of transition: the metal-ordinary insulator transition and the 
mctal-QSH insulator transition. 1 

The cross section of the network is measured in terms of the number of links; a slice of width L contains L links, 
and each link has a spin degree of freedom. In this paper we focus on the case of L even with RBC. (For L odd, 
perfectly conducting channels appear! 13 ' 14 ) 



B. Lyapunov exponents 

We consider a quasi-one-dimensional system. The Lyapunov exponents are estimated using a transfer matrix 
method. The transfer matrix T x relates the current amplitudes ip x on the slice x to those on the next slice, 

ip x+ x = T x i> x . (3) 

Consider the matrix T^T, where T is the transfer matrix product, 

T = T M ---T 2 T 1 . (4) 

The matrix T^T is Hcrmitian and positive definite (because T is invertible) and so its eigenvalues are real and positive. 
They also occur in reciprocal pairs because of current conservation. From each eigenvalue a Lyapunov exponent is 



The direct transition from the ordinary to QSH insulator occurs only for q — > 0. In this limit, the transition point belongs to the same 
universality class as the integer quantum Hall system. 



Critical exponent for the quantum spin Hall transition in Z2 network model 



3 



denned by the limit 



7 = lim , (5) 



where A is an eigenvalue of T'T. In QSH systems, states are doubly degenerate (Kramers degeneracy) and the same 
degeneracy occurs in the Lyapunov exponents. Keeping this in mind, and putting the exponents in decreasing order, 
we number them as 

7l/2 > Tl/2-i > • • • > 72 > 7i > > -71 > -72 > • • • > -Jl/2 • (6) 

Here, it is to be understood that each exponent is doubly degenerate. 

The calculation of the Lyapunov exponent Eq. ([5]) is terminated at finite M when the target Lyapunov exponent 
has converged to within a specified precision. Note that higher Lyapunov exponents converge more quickly than lower 
Lyapunov exponents. 

III. RESULTS 
A. Finite size scaling 

We consider a Z2 network model with strip geometry with RBC in transverse direction and the spin-mixing pa- 
rameter set to q = 0.309. For each F, we assume the scaling formula, 

F(L,p) = F((p - p c )LV», (p - p c )L*"), (7) 

where j/i rr (< 0) is the exponent for the irrelevant correction to scaling. We expand the scaling function around the 
QSH transition point p = p c as 



n d (i) 

:e 

j=l i=l j=0 



T(L,p) = F c + £ Cj (uL 1 '") 3 + E E d ^)( uLl/u y ( vLm "Y > (8) 
where u is the relevant scaling variable, 



and v the irrelevant scaling variable, 



» b 



EM^) - * = 1 > w 



«=EM^ ' 6o = L (10) 

First, we analyze V, which is calculated from the smallest positive Lyapunov exponent (see Fig. [2]). A fit of the 
numerical data to Eq. ([8|) with n a = 3, n& = l, n c = 5, rid(l) = 5, n,d(2) = 2 and taking p c , T c , aj, bj, Cj, dj(i), v, and y- lrr 
as fitting parameters, yielded 

p c = 0.561±0.001, r c =0.1400±0.0004, ^ = 2.64±0.06, j/ lrr = -0.94±0.04. (11) 

The precision of the data for the smallest positive Lyapunov exponents is better than 0.03% for L = 16, 24, 32, 40, 48, 64 
and 0.05% for L = 96, 128, 192, and goodness of fit is 0.2. Nevertheless, the precision of the estimate of the critical 
exponent v is not sufficient to distinguish it from the critical exponent of the quantum Hall transition, v = 2.59±0.01. 15 
Next we analyze T 2 , which is calculated from the 2nd smallest positive Lyapunov exponent (see Fig. [3]). We see 
a clear common crossing point of the curves separating the metallic and topological insulating phases. A fit of the 
numerical data to Eq. ([7J with n a = 2, n(, = l, n c = 4, nrf(l) = 4, n£;(2) = 2 yielded 

p c = 0.562±0.001, r c = 1.429±0.004, ^ = 2.73±0.02, y irr = -0.95±0.02 . (12) 

The precision of the data for the 2nd Lyapunov exponents is better than 0.01%, and goodness of fit is 0.7. This result 
is in good agreement with that obtained by Asada et al*& for SU(2) model with PBC, v — 2.746 ± 0.009, and by 
Obuse et al£ for Z 2 network model with RBC at the metal-ordinary insulator transition, v = 2.88 ± 0.04. It is also 
clearly different from the exponent for the quantum Hall transition. 
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FIG. 2: r of Eq. (JlJ versus tunneling probability for the Z2 network model near the metal-QSH insulator transition. The lines 
are a finite size scaling fit, with different lines corresponding to different cross sections (L = 16, 24, 32, 40, 48, 64, 96, 128, 192). 
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FIG. 3: T2 of Eq. ([2]) versus tunneling probability for the Z2 network model near the metal-QSH insulator transition. The 
lines are a finite size scaling fit, with different lines corresponding to different cross sections (L — 32, 40, 48, 64, 96, f28, 160, 192). 



IV. SUMMARY 



We have reported an estimate the critical exponent for the divergence of the localization length at the metal- 
quantum spin Hall insulator transition. By analyzing the scaling of the 2nd smallest positive Lyapunov exponent, we 
have estimated the critical exponent to be v — 2.73 ± 0.02. Our result shows that this critical exponent is insensitive 
to the topological property of the insulating phase. This is in sharp contrast to the critical conductance distribution, 4 
which is sensitive to the type of transition, i.e., the presence or absence of edge states in the adjacent insulating phase. 

Analysis of the scaling of higher Lyapunov exponents may also be useful in the study of other Anderson transitions 
where conducting edge or surface states occur in the insulating phase, such as, for example, in topological insulators. 
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